Non-Hamiltonian dynamics in optical microcavities resulting from wave-inspired 

corrections to geometric optics 
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We introduce and investigate billiard systems with an adjusted ray dynamics that accounts for 
modifications of the conventional reflection of rays due to universal wave effects. We show that 
even small modifications of the specular reflection law have dramatic consequences on the phase 
space of classical billiards. These include the creation of regions of non-Hamiltonian dynamics, 
the breakdown of symmetries, and changes in the stability and morphology of periodic orbits. 
Focusing on optical microcavities, we show that our adjusted dynamics provides the missing ray 
counterpart to previously observed wave phenomena and we describe how to observe its signatures 
in experiments. Our findings also apply to acoustic and ultrasound waves and are important in 
all situations where wavelengths are comparable to system sizes, an increasingly likely situation 
considering the systematic reduction of the size of electronic and photonic devices. 

PACS numbers: 03.65.-w, 03.65.Sq, 05.45.-a, 40., 91. 
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Quantum-classical and ray-wave correspondence are 
fundamental in achieving a joint understanding of the 
quantum (wave) and classical (ray) worlds. In recent 
years, these ideas have been brought forward especially 
in the field of quantum chaos [l[ of electronic [|J and 
optical systems [||. Besides this fundamental, there is 
a practical aspect: If true, the handy ray picture readily 
provides information about the wave propagation in opti- 
cal cavities, at seismic interfaces, or in acoustics. Indeed, 
ray- wave correspondence was confirmed in numerous sit- 
uations [E 0, E 11, 0, II, > even away from the formal 
classical limit of vanishing wavelength A. In practice, ray- 
wave correspondence is usually taken for granted. How- 
ever, examples of its violation are accumulating recently, 
especially in open systems where A is comparable to the 
system size L or radii of curvature R associated with 
bent interfaces or obstacles. For instance, in optical mi- 
crocavities of spiral shape, "quasiscar" -resonances have 
been computed for which no ray analogue could be iden- 
tified fiol j. Often, Husimi projections [ll[ systematically 
show less symmetry and small mismatches in comparison 
to the classical Poincare surface of section [12, [l3[ . Note 
that discrepancies between ray and wave pictures may 
arise for all types of waves when R, L 2; A. 

In this Letter, we show that these discrepancies arise 
when ignoring deviations from the conventional ray pic- 
ture occurring towards the wave regime. We incorporate 
these deviations in an adjusted reflection law, schemati- 
cally illustrated in Fig. [TJ and we show the dramatic con- 
sequences to the long time dynamics of billiards. Most 
remarkably, the ray dynamics typically becomes non- 
Hamiltonian and, e.g., attractors and repellors replace 
the so-called Kolmogorov-Arnold-Moser (KAM) islands 
of stability. We show how our adjusted ray dynamics ex- 
plains the previous observations of quasiscar-resonances 
and of broken symmetries, and we describe how to test 
our predictions in the laboratory. 

When a light ray hits a glass-air interface, it will split 
into two rays, one reflected and one refracted, unless the 
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FIG. 1: (Color online) Generalized reflection law that includes 
a shift As along the boundary coordinate s (GHs) and an 
increase in the angle of the reflected ray from to 8 + Ad 
(Ff). The arrows denote the beam center. 



angle of incidence Oi (measured to the boundary normal) 
is larger than a critical angle, 9 C = arcsinl/n (where n 
is the refractive index of the material, n — 1.5 for glass) 
and total internal reflection occurs. Reflection and re- 
fraction follow the well-known Snell and Fresnel laws; in 
particular, the reflection angle r will equal the incident 
angle Oi. This is the well-established ray picture of geo- 
metric optics. In reality, however, the incident light ray 
will be a beam, a bundle of rays with angles of incidence 
6*b centered around a mean angle Oi. When Oi « C the 
larger Ob experience total internal reflection, whereas the 
smaller Ob are refracted. Consequently, the mean angle of 
the reflected beam will be larger than the original angle 
of incidence, r > Oi: The ray picture is violated! This 
effect has been termed Fresnel filtering (Ff) [Til ITsl liH ] 
and was recently measured with microwaves |l7l| Re- 
markably, a similar effect was observed in chemical reac- 
tions in excitable media fill ]. 

Another feature inherent to total internal reflection of 
beams, also shown in Fig. [TJ is the Goos-Hanchen shift 
(GHs), discovered in 1947 for light but also well-known 
in acoustics [lS, 20] . Its origin is the interference of rays 
with slightly different angles of incidence that can be in- 
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FIG. 2: (Color online) (a) Spiral billiard with e = 0.1: conventional (dashed line) and adjusted (solid line) ray dynamics 
with constant Ff A9 — 1.737°. (b) Values (4>,A9) for which a period 3 triangular orbit (black line) and a period 5 star orbit 
(red/gray line) exist. The dashed line corresponds to 2a, where a shown in (a) is the small angle between the normal and the 
vector pointing to R = 0, given by sin(a) = e/ ^/e 2 + (1 + e(f>) 2 . 



terpreted as a lateral shift As of the reflected beam along 
the material interface, parametrized by s, the arclength 
of the interface. As both Ff and GHs have the same ori- 
gin (generic beams consisting of a distribution of 0b's) 
they arise together [l|| [l7|, as shown in Fig. [T] but act 
in complementary phase space directions. It remains an 
important open problem to obtain universal analytical 
expressions of the GHs and Ff at curved interfaces. For 
reflections off circular boundaries, the GH and Ff effects 
were computed in [1 61 ] (see Fig.[3]below), where the single 
reflection of a Gaussian beam at a circular air hole scat- 
terer in a glass matrix was carefully investigated. For 
a size parameter nkR = 400 (k — 2ir/\), correspond- 
ing to the semiclassical regime between the pure wave 
and ray limits, Ff and GHs amount to a few degrees. 
They are even larger at smaller nkR and appreciable es- 
pecially near 9 C . In this paper, we adopt the terminol- 
ogy used in optics and refer to the corrections leading 
to non-specular reflection as Ff and GHs, but emphasize 
the potential importance of corresponding effects in the 
reflection of acoustic, seismic [UJ, or even shock [22[ and 
chemical [18| waves. 

Interestingly, the combined effect of Ff and GHs have 
never been incorporated in the ray dynamics of billiards, 
despite their popularity as model systems for example in 
the field of quantum chaos. Here we fill this gap and use 
the results obtained in 

Ref. [l(| for a single reflection to investigate the long- 
term dynamics in billiards. Changes in the phase space 
are expected [HI, not least because of the results 
obtained for soft- wall billiards [HI, UHl • We shall see in 
the following that taking into account Ff and GHs will 
indeed tremendously affect phase space of optical micro- 



cavities that are described with the adjusted ray optics. 
In particular, we will find that the emerging new features 
are in much better agreement with previously obtained 
wave dynamical results - using the adjusted ray optics 
profoundly deepens our understanding of the wave dy- 
namics in open (optical) quantum-chaotic systems. 

We first consider the effect of the adjusted ray dynam- 
ics in spiral billiards defined by R((f>) = 1 + e0/27r, in 
polar coordinates, as shown in Fig. [Ha). We focus on 
rays rotating in the clockwise direction, as the one de- 
noted by the dashed line, because the other rays quickly 
hit the "notch" at 4> — and leave the cavity. For these 
rays, the angular momentum M (being the shortest dis- 
tance from the ray to the center of the spiral) decreases 
with every collision due to the small angle a, described in 
Fig. O This implies that no periodic orbit avoiding the 
notch exists. Ray-wave correspondence suggests that the 
same should be true for (wave) resonances. Therefore, 
the finding of triangular (n = 2) and star-shaped (n = 3) 
resonances in Ref. [lOj which were named "quasiscars" 
came as a true surprise. 

Let us now correct the conventional ray picture by 
adding (constant) Ff and GHs. The Ff increases the 
reflection angle by A9, which counteracts the above- 
mentioned decrease of M. The GHs A(f> increases M, 
however, this effect is roughly e/27t times smaller than 
Ff (for AO A(f>), and will be neglected [23[. For a 
full interval of values AO (around 2a) the conservation 
of M can be fully restored after several collisions, allow- 
ing periodic orbits to exist. The solid curve in Fig. [2ja) 
shows one of such (unstable) periodic orbits, obtained for 
A0 = 1.737°(= 2a for e = 0.1 and = tt). As dis- 
cussed above and shown in Fig. [3] this value is consistent 
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FIG. 3: (Color online) Terms of the adjusted reflection law: GHs As (upper curve) and Ff A9 (lower curve) for a circular 
obstacle with n = 1.5 and kR — 400 (from _lf|). The correction is largest near the critical angle 8 C ; we assume a monotonic 
increase of Ff towards 8 C . Insets: Annular billiard geometry (left) and its conventional phase space (right), obtained through 
a Poincare section (sin# vs. spatial position s) at the inner boundary. 



with the value of Ff expected for the parameters used in 
Ref. OH (6 » 6 C and nkR « 100). Figure[^b) shows that 
changing Ad yields triangular orbits slightly rotated in 
space (an effect also seen in Ref. 10]); in fact, for almost 
all reflection points 4> there is a Ad for which a periodic 
orbit exists. We found equivalent results for star shaped 
orbits and other periodic orbits with small period. The 
authors of Ref. [l(| correctly argued the closeness of di of 
the regular resonances to d c to be crucial. We have now 
clarified that Ff provides the explanation - note that GHs 
alone does not allow the angular momentum M to retain 
its initial value along the orbit - and that the regular 
trajectories seen are thus real scars [30] of a spiral bil- 
liard with adjusted ray dynamics. The above results lead 
to the following predictions for the experiments that are 
being performed: A6(k) decreases as kr increases, and 
so quasiscars are expected to appear at different loca- 
tions <f> until Ff becomes too small and no quasiscars are 
predicted to appear, in agreement with the semiclassical 
limit. We predict also that depending on n other orbits 
will be stabilized, e.g., the square orbit for n w y2. 

Let us now investigate in more detail the joint effect 
of Ff and GHs on the phase-space dynamics (s,sin6>), 
that is most often used to establish the ray-wave corre- 
spondence. We adjust the conventional reflection law by 
replacing 



S{ + As(si,9i) 
Oi + Adisi^i) 



(1) 



The dynamics in a billiard is described by the billiard 
transformation T that maps a collision (si,sin(?i) to the 
next one. Considering ((T|), the absolute value of the Ja- 
cobian of T is given by 
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To understand the new billiard dynamics is crucial to 
know whether J deviates from 1. The phase-space area 
cos ddsdd is preserved if and only if J = 1, which, 
in two dimensions, implies Hamiltonian dynamics (ev- 
idently, J = 1 for As = Ad = 0). For J ^ 1, the 
map does not preserve the area cos ddsdd and typically 
non-Hamiltonian dynamics arises. Equation (|2|) indicates 
that J sensitively depends on the system-specific depen- 
dence of A8 and As on 8 and s, allowing, in principle, J 
to take arbitrary values. 

We are interested on how the modification ((T|) affects 
the dynamics of a billiard for which the conventional re- 
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FIG. 4: (Color online) Phase space of the annular billiard [2||] with (a) conventional reflection law (magnification of right inset 
of Fig. [3j, and including (b) GHs, (c) Ff, and (d) GHs and Ff, as described in Fig. [3j Note that GHs modifies KAM islands 
in (b), but preserves the Hamiltonian character that is lost when Ff is taken into account in (c) and (d), where attractors (a) 
and repellors (o) are formed. Insets show the corresponding orbits in the real space. 



flection law produces a generic phase space, i.e., coex- 
istence of chaotic and regular dynamics. A long stud- 
ied [28| and physically realistic example is the annular 
billiard shown in the left inset of Fig. [3] A glass disk with 
conventional reflection law at the outer boundary (real- 
izable by a mirror coating) and an eccentrically placed 
inner disk scatterer (of radius r and distance S between 
the disk centers). Ff and GHs act thus only at the in- 
ner boundary, where the Poincare section is defined. For 
GHs and Ff we use the values reported in Ref. |l6[, cf. 
Fig. [3] We are not interested in the refractive intensity 
loss at the inner boundary, that could be included via 
Fresnel's law [29[. Rather, our focus is on the dynamics 
of rays that remain in the annulus. 

In this system, the curvature is constant, R = r, and 
Eq. ||2J| simplifies to 



J = 



cos 



A9) 



dA9" 



cos(0) V 1+ de ) ■ (3) 

Because A9 is a small correction, mainly the last term 
dA9/d9 determines J. Three regions, with J > 1, J < 1, 



and J = 1, are marked in Fig. [3J Next we compare 
the phase space of billiards with conventional dynamics 
and those with GHs and Ff corrections. We have chosen 
parameters r — 0.2, S — 0.3 in the annular billiard such 
that at least one KAM island is seen in each of the three 
regions, named II, 12, and 13 and indicated in the right 
inset of Fig. [3l The boxed region is enlarged in Fig. 
Hk, where the stable periodic orbits in the center of the 
islands are shown as insets. The following panels (b)-(d) 
show that the effect of GHs and Ff can destroy, create, 
and modify these KAM islands. Now, we explain these 
observations. 

The change in the structure of the phase space is re- 
lated to the value of J, cf. Eq. {3J and Fig. g] While 
elliptic islands contained in the area-preserving region 
J = 1 {9 3> 9 C , II) survive when GHs and Ff are switched 
on, those contained in the region J < 1 (J > 1), where 
the dynamics is area contracting (expanding), become 
attractors (repellors) once the Ff correction is taken into 
account. The area contracting (J < 1) region is placed 
around the critical angle 8 w 8 C and affects 12 that gives 
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rise to a periodic attractor in Fig. [4jc)-(d). 

Note that although this (effective) dissipation is re- 
lated to the openness of the system, via the partial re- 
fractive escape of light that gives rise to the Ff effect, 
the refractive escape itself is not its origin. This becomes 
evident as also J > 1 is realized for subcritical incidence 
9 < 9 C , and 13 originates a repellors in Fig.[4jc)-(d). This 
dissipative behavior is related to the openness of the sys- 
tem via the partial refractive escape of light that gives 
rise to the Ff effect. Note, however, that the possibility 
of refractive escape and the Ff effect may also generate 
expanding dynamics (J > 1), i.e., the first-sight impres- 
sion that refractive escape can only lead to contracting 
dynamics because of the loss of energy through trans- 
mitted rays does not apply. (Note also that we are not 
interested in the refractive intensity losses in the present 
study.) Whether the Jacobian J takes a value larger or 
smaller than 1 depends, in the case of the annular bil- 
liard, cf. Eq. ([3]), on dA6/86, i.e., on the dependence of 
the Ff correction A9(9) on the angle of incidence 9 [see 
Eq. ||2J| for the general case]. Indeed, for subcritical in- 
cidence 9 < 6> c we find dA9/d9 > and Eq. (3) leads to 
J > 1. Accordingly, the KAM island 13, present in the re- 
gion of subcritical incidence in Fig. 4(a), transforms into 
a repellor when Ff is considered in Fig. 4(c)- (d). The 
action of GHs alone, cf. Eq. ([3]) and Fig. [4p, preserves 
the Hamiltonian dynamics (J — 1). In this case, mod- 
ifications of the phase space similar to those obtained 
by changing a control parameter occur, e.g., in Fig. 2J 
region 12 is absent and small islands are created in the 
secondary chain of II and 13 (hardly visible) . Most excit- 
ing and dramatic in its effect is the Ff contribution that 
renders islands into attractors (close to 9 C ) and repellors, 
i.e., signatures of a non- Hamiltonian dynamics (J =/= 1) 
are seen in a model system often used in quantum and 
wave chaos. 

Another important effect of GHs and Ff is the break- 
ing of symmetries. The conventional phase space of bil- 
liards with one symmetry axis, like the annular billiard, 
possesses the time reversal symmetry I t : (s,sin#) — > 
(s, — sin#) as well as a spatial symmetry I s : (s, sin#) — > 
(2ir — s, — sin (9). In the example above, Ff violates the 
it-symmetry. The simplest way to see this is to consider 
a beam with (average) critical incidence. The subcritical 
part of the beam will be refracted, and the remainder will 
be reflected to supercritical angles, violating the rules of 
specular reflection. If this reflected part of the beam is 
now considered as the (time-reversal partner of the origi- 
nal) incident beam, it becomes immediately evident that 
time-reversal symmetry and the principle of ray-path re- 
versibility are lost as all rays of this beam will be totally 
reflected and specular reflection is obeyed. As a conse- 
quence of the violation of time-reversal symmetry due to 
Fresnel filtering, the I t symmetric orbit in the center of 
12 shown in Fig. |4ja) becomes Ff- and GHs-distorted, 
see the left inset in Fig. 0Jd). Moreover, for the oppo- 
site sense of traveling (sin# < 0), the distortion will be 
in the opposite direction. This shows that the position 



of the reflection point at the outer boundary depends on 
the sign of 9, and that /(-symmetry is broken. In wave 
calculations for optical microcavities (where Ff and GHs 
are naturally included), indications of /(-breaking in the 
Husimi projection were noticed [HI, EH, but their ori- 
gin remained unclear. The adjusted reflection law that 
we propose here obviously accounts for this symmetry 
loss. 

The consequences of an adjusted, non-specular reflec- 
tion law are accessible to experiments using, e.g., state of 
the art microwave resonators with the geometry adjusted 
such that periodic orbits (KAM islands) exist close to 
the critical line. For example, the annular geometry sug- 
gested here can be realized as a teflon disk with an air 
hole and an outer metallic coating. One approach is the 
analysis of the length spectra, taken as Fourier transform 
of the resonances in a finite interval of wavenumbers k. 
For different k intervals, the resonant wave patterns ad- 
just to the new fci?-dependent Ff and GHs. Therefore we 
expect changes in the peak position (length of the ad- 
justed orbit) and peak height (stability of the adjusted 
orbit or size of KAM island) . Another possibility is to di- 
rectly observe the /c-dependent morphology of resonances 
by directly measuring their wave patterns. 

In principle, our results also apply to ballistic elec- 
tronic mesoscopic systems such as quantum dot billiards 
[32| since the Fermi wavelength of the electrons is lower, 
but comparable to the system size. In addition, the con- 
finement of the electrons in the two-dimensional electron 
gas layer is by gate electrodes and will never be perfect. 
In this sense the confining potential has soft walls, and 
corrections to the hard-wall ray dynamics are expected. 
We point out that most works in electronic mesoscopic 
systems focus on transport properties. Such an openness 
of the system (by attaching leads and applying bias volt- 
ages) is, however, very different from the openness that 
we consider here. 

In conclusion, we have investigated billiards with an 
adjusted ray dynamics that incorporates wave correc- 
tions to the law of specular reflection. We have shown 
that the conventional phase space experiences dramatic 
changes, including the creation of periodic orbits, the 
breakdown of symmetries, and the creation of regions of 
non-Hamiltonian dynamics where attractors and repel- 
lors replace KAM islands. We have presented strong 
evidences that these results explain previous observa- 
tions of 'quasi-scared' modes in spiral microcavities and 
of broken time reversible symmetry. We expect that 
other violations of the ray-wave correspondence can 
be explained similarly, e.g., the observation of dispro- 
portionally man y r egular, scarred,and scar-like reso- 
nances 0, i, [Mill M, EE [M M, El might be related 
to the creation of periodic attractors, and discrepancies 
in the far-field pattern [l2|, HH might be explained by 
the effect of Ff and GHs on invariant sets and their un- 
stable manifolds. Although we focused on the case of 
optical microcavities, where Ff and GHs are well estab- 
lished effects, our main results apply likewise to acous- 



6 



tic, seismic, and ultrasound waves whenever deviations 
from the specular reflection become relevant, e.g., when 
the wave length becomes comparable to the system size. 
We want to emphasize the importance of previous un- 
explained disagreements between ray and wave results, 
including those potentially regarded as numerical or ex- 
perimental errors, because they might be explained by 
the adapted ray dynamics introduced here. 
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